1 



THE STRUCTURE OF THE PHOTON* 



STEFAN SOLDNER-REMBOLD 
Albert- Ludwigs-Universitat Freiburg, Hermann- Herder-Str. 3 
D-79104 Freiburg, Germany 
E-mail: soldner@ruhpb.physik.uni-freiburg.de 

The structure of the photon is studied in high energy photon-proton interactions at HERA and 
photon-photon interactions at LEP. The status of these measurements is reviewed. 
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1 Introduction 



The photon is one of the fundamental gauge bosons of the Standard Model without self- 
couplings and without intrinsic structure. However, at high energies photon-hadron and 
photon-photon interactions are dominated by quantum fluctuations of the photons into 
fermion-antifermion pairs and into vector mesons which have the same spin-parity ( J PC = 
£<\ • 1 ) as the photon. This is called photon structure. Electron-positron collisions at LEP and 

positron-proton collisions at HERA are an ideal laboratory for studying photon structure 
in the interactions of quasi-real and virtual photons, testing predictions of both Quantum 



' Electrodynamics (QED) and Quantum Chromodynamics (QCD). 

t*"" ' 2 Electron-photon scattering 

OV 

If one of the scattered electrons in e + e~ collisions is detected (tagged), the process e + e~ — > 
e + e~ + hadrons (Fig. p]) can be regarded as deep-inelastic scattering of an electron^] on a 
Qh. quasi-real photon which has been radiated by the other electron beam. The cross-section is 

written as 



d Ce-y— > e +hadrons 27T(J 



l + (l-y) 2 )F2(x,Q 2 )-y 2 FZ(x,Q 2 ) , (1) 



|Xj ■ dxdQ 2 ~ xQ 4 

. . , where a is the fine structure constant and 

Q 2 = -q 2 = -{k-k') 2 

is the negative four-momentum squared of the virtual photon 7* and 

Q 2 Q 2 

X 2p-q Q 2 + W 2 + P 2 
p-q 



V 



p ■ k 



'Invited talk given at the XVIII International Symposium on Lepton Photon Interactions, Hamburg, Ger- 
many, July 28-August 1, 1997 

"In this paper positrons are also referred to as electrons 



FREIBURG-EHEP-97-19 



2 



are the usual dimensionless variables of deep-inelastic scattering. W 2 = (q + p) 2 is the 
squared invariant mass of the hadronic final state. The negative four-momentum squared, 
P 2 = —p 2 , of the quasi-real target photon is approximately zero and therefore usually 
neglected. In leading order (LO) the photon structure function F^x^Q 2 ) is related to the 
sum over the quark densities of the photon weighted by the quark charge e q 



F2(x,Q 2 ) 



2*£ 



^q/q/7 i^i Q ) 



with fq/j(x, Q 2 ) being the probability to find a quark flavour q with the momentum fraction 
x (sometimes denoted by x 7 ) in the photon. For measuring F^{x, Q 2 ) the values of Q 2 and y 

can be reconstructed from the energy, £tag, and the 
angle, # ta g, of the tagged electron and the beam 
energy -Ebeam via the relations 

Q 2 2£ , beam £ , ta g(l — COS^tag) 




w 2 



y 



1 



E, 



tag 



COS 



Ttag 

2 



anti-tag (-p 2 =P small) 



-^beam 

In order to identify an electron in the detector, the 
tag energy E tSbg has to be large, i.e. only low values 
of y are accessible (y 2 <C 1). The contribution of the 
cross-section term proportional to the longitudinal 
structure function F? is therefore negligible. 



Figure 1: Deep-inelastic e*y scattering: k(k') 
denotes the four-momentum of the incom- 
ing (scattered) electron and q(p) the four- 
momentum of the virtual (quasi-real) photon. 

The reconstruction of x, however, relies heavily 
on the measurement of the invariant mass W from the energies E^ and momenta ph of the 
final state hadrons h: 

W 2 = 

Unfolding of the x dependence of the structure function therefore requires that the hadronic 
final state in events is well measured and well simulated by the Monte Carlo models. 





2. 1 Hadronic energy flows 

OPAlI, ALEPrl! and DELPHli have studied the hadronic energy flow per event, 1/N ■ 
dE/drj, as a function of the pseudorapidity 77 = — In tan 8/2, where the sign of 77 is chosen 
in such way that the tag electron is always at negative rj. OPAL has measured the energy 
flow at medium Q 2 ((Q 2 ) = 13 GeV 2 ) and at high Q 2 ((Q 2 ) = 135 GeV 2 )0. In Fig. the 
energy flows are compared to the two QCD based Monte Carlo generators HERWIGcl and 
PYTHIaI. The data distributions have been corrected for detector effects. The generator 
F2GEN is used to model the unphysically extreme case of a two-quark state in the 7*7 
centre-of-mass system with an angular distribution as in lepton pair production from two 
real photons ("pointlike"). The "perimiss" sample is a physics motivated mixture of point- 
like and peripheral interactions, where peripheral means that the transverse momentum of 




Figure 2: The energy flow per event as a function of the pseudorapidity 77 in different Q 2 ranges compared 
to the HERWIG, PYTHIA and F2GEN Monte Carlo models. The data have been corrected for detector 
effects. 



the outgoing quarks is given by an exponential distribution as if all the photons interacted 
as pre-existing hadrons. 

Significant discrepancies exist between the data and all of the Monte Carlo models. The 
agreement improves at higher Q 2 . Since x and Q 2 are correlated, the discrepancies at low 
Q 2 are observed also at low x. These discrepancies between the data and the Monte Carlo 
model for the hadronic_final state are the dominant source of systematic uncertainty in the 
unfolding of 



x 



Q 2 )l 



Tuning of the Monte Carlo to improve the mod- 
elling of the hadronic final state is complicated and 
must be done with care in order to avoid a bias of 
the result of the unfolding towards the parametri- 
sation of the parton distribution functions used in 
the tuned Monte Carlo. ALEPH0 has measured the 
energy flow in tagged event for (Q 2 ) = 14.2 GeV 2 . 
The distributions have not been corrected for de- 
tector effects. The energy flow shown in Fig. || 
is compared to the HERWIG generator @ and a 
Monte Carlo which consists of a mixture of Quark 
Parton Model (QPM) and Vector Meson Domi- 
nance (VDM) similar to the F2GEN generator. In 
addition, a modified version of HERWIG ("HER- 
WIG+power law pr") was used. The modifica- 
tion is based on studies of the photon remnant by 
ZEUsi. In standard HERWIG a Gaussian distribu- 
tion is used to describe the limitation of the trans- 
verse momentum of the outgoing partons with re- 
spect to the initial target photon. In the modified 



3.5 



-f ALEPH data 
HERWIG 5.9 

— HERWIG 5.9 + power law pt,_ 

— QPM + VMD 




Figure 3: The energy flow per event as a 
function of 77, compared to a QPM+VDM 
model, the standard and the tuned version 
of HERWIG. 
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HERWIG the Gaussian is replaced by a power law spectrum. The agreement with the data 
improves. A similar study was performed earlier in Ref. 7 using OPAL data. It is expected 
that such improvements of the Monte Carlo models will significantly reduce the systematic 
error of the structure function measurements for hadronic events. 



2.2 The photon structure function at high Q 2 

Even though the concept of the photon structure function i 7 ^ has been developed in analogy 
to the formalism of the nucleon structure functions , there are important differences: 
F^{x, Q 2 ) increases with Q 2 for all x and this positive scaling violation is expected already 



within the parton model. Furthermore, is large for high x, whereas F<2 decreases at 
large x. These differences are due to the additional perturbative 7 —* qq splitting which 
does not exist for the nucleon. 

For large x and asymptotically large Q 2 the value of can therefore be calculated 
from perturbative QCDtJ. The next-to- leading order (NLO) result IHI can be written as 



El 

a 



a(x) 



+ b(x) 



(2) 



where a(x) and b(x) are calculable functions which diverge for x — > and a s As the strong 
coupling constant. The first term corresponds to the LO result by Witten Q. The mea- 
surement of i*2 could be a direct measurement of Aqcd if it were not for the large non- 
perturbative contributions due to bound states. These contributions are large at all exper- 
imentally accessible values of Q 2 . 



> 
o 



gj Vs=91 GeV (LEP1) 
; Vs=180 GeV (LEP2) 






Q 2 ='W,?„*/(l-x) 














Q 2 =<#/(l-x) - 
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+ OPAL (0.1 < x < 0.6) 
if AMY (0.3 < x < 0.8) 
t JADE (0.1 < x < 1.0) 
\ DELPHI prel. (0.3 < x < 0.8) 
if TPC (0.3 < x < 0.6) 

— GRV LO (0.1 < x < 0.6) 
GRV LO (0.2 < x < 0.9) 

— - GRV LO (0.3 < x < 0.8) 

SaSID (0.1 < x < 0.6) 
HO (0.1 < x < 0.6) 

— ASYM(0.1<x<0.6) 



t TOPAZ (0.3 < x < 0.8) 
if ALEPH prel. (0.3 < x < 0.8) 




Q- [GeV ] 



Figure 4: a) Kinematical coverage of the (Q 2 ,x) plane at LEP1 and LEP2. b) The photon structure 
function Fjj! /a as a function of Q 2 . 

The photon structure function F%(x, Q 2 ) can be measured at LEP in the range x > 10~ 3 
and 1 < Q 2 < 10 3 GeV 2 which makes it possible to study the QCD evolution of in a 
wide range of x and Q 2 (Fig. ^|a). 
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The evolution of F^ with InQ 2 is shown in Fig. |4]b using the currently available K7 
measurements for 4 active flavours. The data are compared to the LO GRVcJ and the 
SaS-lD0 parametrisations, and to a higher order (HO) prediction based on the NLO GRV 
parametrisation for light quarks and on the NLO charm contribution calculated in Ref. 12. 
The data are measured in different x ranges. The comparison of the LO GRV curves for 
these x ranges shows that for Q 2 > 100 GeV 2 significant differences are expected. An 
augmented asymptotic prediction for F^ is also shown. The contribution to F^ from the 
three light flavours is approximated by Witten's leading order asymptotic form El. This 
has been augmented by adding a charm contribution evaluated from the Bethe-Heitler 
formulaEH, and an estimate of the hadronic part of F%, which essentially corresponds to the 
hadronic part of the LO GRV parametrisation. In the region of medium x values studied 
here this asymptotic prediction in general lies higher than the GRV and SaS predictions but 
it is still in agreement with the data. The importance of the hadronic contribution to F% 
decreases with increasing x and Q 2 , and it accounts for only 15 % of at Q 2 = 59 GeV 2 
and x = 0.5. 

As predicted by QCD the evolution of leads to a logarithmic rise with Q 2 , but 
theoretical and experimental uncertainties are currently too large for a precision test of 
perturbative QCD. 



GRV-NLO - -- SaS-lD (LO) 



sr 0.5 

0.25 


0.5 

0.25 


1 

0.5 


1 

0.5 




PLUTO 9.2 
OALEPH prl 8.' 
▼ OPAL 9.0 



TPC/2y 
Q 2 =0.24 GeV 2 



_1_ 



TPC/2y 2.8 



_L_ 



DELPHI 12 
(A DELPHI prl 1 



OPAL 135 




TPC/2y 0.38 



OPAL 3.76 



PLUTO 4.3 



▼ OPAL 14.7 
J&OPAL 14.5 



_L_ 



ALEPH prl 279 AMY 390 



TPC/2y 0.71 




TPC/2Y 1.3 



>TPC/2y 5.1 
I TOPAZ 5.1 



DELPHI prl 6 
• AMY 6.8 



TOPAZ 16 




OPAL 1.86 



ALEPH prl 19. 



PLUTO 2.4 



30PAL 7.5 



DELPHI prl 2|2 
OTASSO 23 
A JADE 24 



1 

X 



Figure 5: Measurements of the photon structure function in bins of x and Q 
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2.3 The photon structure function at low x 

All currently available measurements^ of the photon structure function are shown in Fig. || 
The data are compared to the NLO GRV0 parametrisation and the LO SaS-lD0 parametri- 
sation. If the photon is purely hadron-like at low x, a rise of the photon structure function is 
expected in the low x region for not too small Q 2 , similar to the rise of the proton structure 
function observed at HERA. Only with the high statistics and high energy LEP2 data will 
it be possible to access regions in x and Q 2 where the rise of K7 could be observed. An 
interesting new low x measurement of is presented by OPAL in the x and Q 2 ranges 
2.5 x 10~ 3 < x < 0.2 and 1.1 < Q 2 < 6.6 GeV 2 . The measurement is consistent with 
a possible rise within large systematic errors. It should be noted that the OPAL points 
are significantly higher than the previous measurement by TPC/27 in a similar kinematic 
range. 

In the kinematic region of the LEP measure- 
ments, the effect of the P 2 evolution can be- 
come quite important. The ALEPH and OPAL 
F2{x,Q 2 ) measurements have all not been cor- 
rected for this effect, i.e. they are actually mea- 
surements of FiJixjQ 2 ,P 2 )- OPAL estimates 
that the effect of the non-zero virtuality P 2 could 
increase their low x measurement of F% by about 
10 % for P 2 = 0. This estimate is based on the 
P 2 dependent parametrisation of the narton dis- 
tributions by Schuler and Sjostrandlij. 

Since there is little experimental information 
about the low x region, Gurvich, Abramowicz 
and Levy (GAL) have tried to use the F% data 
at low x in order to constrain their new F% 
parametrisation 111. Gribov factorisation is as- 



Q' = 4.3GeV' 


Q' = 5.1 GeV' 


\ Q' = 7.2GeV' 


\ Q' = 9.2GeV' 




3 = 14.7GeV ! 


Q 1 = 1 6GeV ! 


r jP \ 


Q' = 45GcV ; 


Q 2 = 77GcV 


Q' = IDOGcV 1 









Figure 6: The GAL parametrisation (continu- 
ous line) compared to the F% data and the scaled 
F£ data using Gribov factorisation. The dashed 
line is the LO GRV parametrisation and the dot- 
ted line the SaS parametrisation. 
photon induced cross-sections. Furthermore this factorisation assumption should also hold 
for virtual photons at low x. The proton structure function i 7 ^ and the photon structure 
function Fi] are then related by the total 7p and pp cross-sections: 



sumed, based on the idea that at high enough 
centre-of-mass energies all total hadronic cross- 
sections are dominated by an universal Pomeron 
trajectory which allows to relate hadron and 



x,Q 2 



F»(x,Q' 



a 



TP 



(W 2 



PP (W 2 ) ' 

The total cross-section are obtained from a Donnachie-Landshoff parametrisation of the 
data (see Sect. ||). This procedure is used to obtain pseudo-data for at low x. The 
evolution of the parton densities is done using a LO Altarelli-Parisi evolution with a starting 
scale Qq = 4 GeV 2 . The predicted rise is much lower at low x than predicted by the LO 
GRV parametrisation which is also shown in Fig. |(| 
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2.4 Azimuthal correlations 

Only the structure function i 7 ^ has so far been determined directly from measurements of 
double-differential cross-sections for e7 events with hadronic or leptonic final states. It has 
been pointed out 113 that azimuthal correlations in the final-state particles from two-photon 
collisions are sensitive to additional structure functions. AzimuthaLcorrelations can thus 
supplement the direct measurement of structure functions. ALEPHta, L3e3 and OPAL EH 
have measured azimuthal correlations using single-tag e + e~ — * e + e _ /i + /i~ events. 

For single-tag events two independent angles can be defined in the 77* centre-of-mass 
system assuming that the target photon direction is parallel to the beam axis: The azimuthal 
angle x is the angle between the planes defined by the 77* axis and the two-body final state 
and the 77* axis and the tagged electron. The variable r] = cos 6* is defined by the angle 
9* between the \i~ and the 77* axis. 

Neglecting the longitudinal component of the target photon, the cross-section can be 
written as (F 2 7 = 2xF^ + Fl): 



da(e7 — > e// + /U 



2vra 2 / 1 + (1 - yf 



dxdydr/dx/27r Q 2 



xy 



2xF} + Fl 



F 1 

F\ cos x + ~y cos 2x 



(3) 



The conventional structure functions are recovered by integration over r] and x : F? 



f-i Jo 



2?r d\dri p7 / 



2?r 



2, A, B). The structure functions Ft and Fl are related to the scatter- 



ing of transverse and longitudinally polarized virtual photons on a transverse target photon, 
respectively. Fj± is related to the interference terms between longitudinal and transverse 
photons and Fg to the interference term between purely transverse photons. The longi- 
tudinal structure function F^ has been shown to be equal to the structure function Fg in 
leading order and for massless muons, although coming from different helicity states of the 
photons. 



Results for 1994 Data 



Results for 1994 Data 




(b) 



S ALEPH PRELIMINARY 



1/2 F„/F- 




Figure 7: The ratio of structure functions, FJ^/F^ and 1/2 ■ Fg/F^ , for the process e + e — > e + e fl . 
The lines show the QED expectation. The ALEPH measurement is not corrected for the cut \rj\ < 0.8. The 
Q 2 ranges are (Q 2 ) = 8.8 GeV 2 (ALEPH) and 1.4 < Q 2 < 7.6 GeV 2 (L3). 
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The variation of Fj^ and Fg with x is in general consistent with QED (FigjT] and Ref. 21). 
The measured values are significantly different from zero. Apart from being an interesting 
test of QED, these results are especially important as a first step towards measuring the 
additional structure functions for hadronic events using azimuthal correlations. Such a 
measurement will be much more difficult due to the problems related to the jet finding in 
hadronic events. 



3 Jet production in 77 and 7p scattering 

If the virtuality Q 2 of the probing photon in ep scattering or the virtualities Q 2 and P 2 
of both interacting photons in e + e~ scattering are approximately zero, the photons can 
be considered to be quasi-real. HERA can then be used as photon-proton collider with 
7P centre-of-mass energies in the approximate range 40 < W 1P < 300 GeV and LEP2 as 
a photon-photon collider with 77 centre-of-mass energies in the approximate range 10 < 
W 77 < 120 GeV. Anti-tagged 7p events at HERA and anti-tagged 77 events at LEP have 
a median Q 2 of about 10 -4 GeV 2 for Q 2 < 4 GeV 2 . The anti-tagging condition is fulfilled 
if no scattered electron was found in the main calorimeters. Events with Q 2 < 10~ 2 GeV 2 
are tagged using small angle calorimeters. 

In LO different event classes can be defined in 77 and 7p interactions. The photons 
can either interact as bare photons ("direct") or as hadronic fluctuation ("resolved"). The 
events are classified by the photon interaction which leads to the nomenclature shown in 
Fig. ||. In NLO the separation into event classes seizes to be unique, since it depends on 
the factorisation scale for the photon, and only the sum of the cross-section for the different 
event classes is well defined. 

Within the LO picture we expect two hard jets, i.e. with large transverse energy Ej, , in 
the final state in addition to possible photon and/or proton remnants. These jets are related 
to the underlying parton dynamics and can therefore be used to constrain the structure of 
the colliding particles, photons and protons. Direct and resolved events can be separated 
by using the fraction x^ of the photon's momentum participating in the hard interaction. 
The direct events are expected to have x~ = 1. For 7p events a separation of the event 
classes has been achieved experimentally by defining the variable il 

obs _ Z^jets-^T e 

2yE e ' 

with E e = 27.5 GeV being the electron beam energy. The sum runs over the two jets with 
highest Ej 1 and Ey = yE e is the energy of the initial state photon. For two-jet final states, 
this observable is equivalent to the LO definition of x 1 . Fig. ^a shows the uncorrected 
x° bs distribution measured by ZEUS in comparison to the LO MC models PYTHIA and 
HERWIG. The peak at high x° bs due to the direct events is observed in the data. However, 
it is smeared out due to hadronisation, QCD radiation and detector effects. 

The measurement of jet cross-sections allows comparisons with perturbative QCD calcu- 
lations which are based on the matrix elements M(cos8*) for all possible parton scattering 
processes. The various matrix elements depend on the parton scattering angle 9* in the 
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Figure 8: Examples for leading order processes in 77 and 7p interactions. 



parton-parton centre-of-mass frame. A ZEUS measurement of the different angular distribu- 
tions is shown in Fig. ^o. The cos \ 9*\ distribution is here given in the jet-iet centre-of-mass 
system for resolved (x° bs < 0.75) and direct dijet events (x° bs > 0.75) cl. The dijet in- 
variant mass of these events is larger than 23 GeV. The dominant LO diagrams for direct 
7P processes involve quark (fermion) exchange whereas those in resolved processes involve 
gluon (boson) exchange. This leads to different angular distributions which agreewell both 
with the LO QCD curves and with the NLO calculations by Harris and Owens Ej. Similar 
results for 77 scattering were reported by OPAL El. 
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Figure 9: a) Uncorrected a;° bs distribution for 7p dijet events, b) dijet angular distribution dcr/d| cos#*| in 
the jet-jet centre-of-mass frame. 



As we have seen, jet production is sensitive to the parton content of the photon. In the 
kinematic range covered by e + e _ experiments the measurements are mainly probing the 
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quark content of the photon, since the gluon distribution only enters in higher orderf]. This 
is different in 77 and 7p jet production where a large fraction of the events are due to gluon 
initiated processes. Measurements of jet production cross-sections can therefore be used to 
constrain the relatively unknown gluon distribution in the photon. 

Different groups have followed different philosophies for extracting information about 
the parton content of the photon from jet cross-sections. In the first approach hadronic jet 
cross-sections are measured and compared to calculations which use different parametrisa- 
tions of the photon's parton densities as input. In the second approach, which is mainly 
followed by HI, LO parton densities are extracted from the measurements. 

For the comparison of theory and experiment it is very important to define suitable 
jet algorithms. Measurements of jet cross-section in 77 and 7p interactions are usually 
made with cone jet finding algorithms c3, since resolved photon interactions are similar to 
hadron-hadron interactions. A problem appears for overlapping jets which can be merged in 
iterative cone algorithms^ This is not described by a NLO calculation which only contains 
three final state partonsc3. 

Other effects like hadronisation or the underlying event in resolved events are also not 
contained in the theory. Underlying event is a loose term for additional activity which, for 
example, can be caused by soft or hard interactions of the photon and proton remnants of 
the same event. The influence of these effects has to be studied in detail before statements 
about the parton distributions can be made, since the additional energy in the event directly 
affects the jet cross-sections. The Monte Carlo models treat this effect by introducing 
multiple parton interactions (MI). 



3.1 Jet shapes 

The internal structure of jets was studied by ZEUS0 by measuring jet shapes in photopro- 
duction. The jet shape is defined as the average fraction of Ej? that lies inside an inner 
cone of radius r concentric with the jet defining cone: 



where E^(r) is the transverse energy within the inner cone of radius r and iVj e t s the total 
number of jets. The 7^ et dependence of the jet shapes is compared at a fixed value for 
r, ip(r = 0.5), to the PYTHIA prediction with and without multiple interactions (MI) in 



Fig. 10a. The larger colour charge of gluon jets leads to increased QCD radiation which 
broadens the jets, an effect which is well known from jet production in e + e~ scattering. In 
the forward direction (positive rf et ) the fraction of gluon jets in the MC increases, but also 
the influence of MI increases. Both effects are expected to lead to the broadening of the 
jets at high rf et which is observed in the data. In Fig. |l0| b the jet shapes tp(r){Qi different 
ranges are also compared to fixed-order perturbative QCD calculations tfj'Eil. Jets get 
narrower with increasing Ej, . A new parameter R sep was introduced in the calculation to 

^Contrary to the proton case, there exists no simple isum rule for the gluon distribution in the photon. A 
sum rule has been derived by Frankfurt and GurvichCJ and SchulercJ. 
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Figure 10: a) Jet shape ip(r) at a fixed value r = 0.5 as a function of ry jct for jets with > 14 GeV; b) 
jet shapes ip(r) as a function of the radius r for jets with — 1 < rf ct < 2 in different regions. 



mimic the effects of overlapping and merging of jets mentioned in the previous section. Two 
partons are not merged into a single jet if their separation in the ncp plane is more than 
i? sep . This problem only exists for the iterative cone algorithms used here, whereas the fop 
clustering algorithm El requires no additional parameter R sep . The choice i? sep = 1.4 gives 



a reasonable description of the data, but at low E 3 ^ 1 a higher i? sep is needed. This increase 
can be attributed to hadronisation effects and a possible underlying event. 



at low Ej?" 



3.2 Inclusive jet cross-sections and NLO calculations 

NLO jet cross-sections for -yyll'H and 7p interactions US'!! have been calculated by many 
authors. To calculate jet cross-sections in perturbative QCD a hard scale is required in the 
event which is usually the transverse momentum px of the final-state partons (or the jet). 
For the calculation it is assumed that the concept of factorisation can be applied. The LO 
jet cross-section is written as convolution of the parton density /j/ 7 of the photon and, in 
the case of 7p scattering, the parton density fj/ p of the proton with the LO matrix elements 
M{j for the scattering of two partons % and j: 



O-0"ep(ee) v 

dxidxod cos 9* 4r* 



f ih (x up l) f j/PM (x 2 ,p* T ) (cosnf 

xx x 2 



The variables xi,x 2 are either the momentum fractions x~, x p of the partons in the photon 
and the proton (for 7p scattering) or the momentum fractions , x~ of the partons in the 
two interaction photons (for 77 scattering). In addition the photon flux from the electrons 
is taken into account using the Weiszacker- Williams Approximation. 

The NLO correction term depends on the two factorisation scales, the renormalisation 
scale and the cone size R. NLO parton distributions of the photon should be used for 
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a consistent NLO calculation. Several such NLO parametrisations are available GRV0, 
GS96!1 and AFGll. 
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Figure 11: The measured inclusive e + p jet cross-section as a function of rj ,et for jets with > 14 GeV 
and E 1 * > 21 GeV using a cone size R = 1. The kinematic regime is defined by Q 2 < 4 GeV 2 and 
134 < W < 277 GeV. The bands represent the uncertainty due to the calorimeter energy scale. The curves 
are the NLO calculations. 



has compared measured inclusive single jet cross-sections in different ranges of 
and rf et to theNLO calculation by Klasen and Kramer S (Fig. [n|). The NLO GRV0 
and the the GS96E3 parametrisations are used for the photon. The discrepancies are large 



jet 
T ■ 



for smaller Ej, in the forward region rf ct > 0.5. This discrepancy disappears at higher E; 
The agreement also improves if a smaller cone size is used (R = 0.7) which corresponds 
to an effective increase of the Ej, threshold compared to R = 1. In the region of large 
discrepancy non-perturbative contributions from the underlying event are expected to be 
large. 

The NLO calculations are given for two values of i? S ep (-^sep = R, 2i2) which indicates 
part of the theoretical uncertainty. The differences between NLO GRV and GS96 are of 
similar magnitude as the current theoretical and experimental uncertainties. In addition 
ZEUS has also measured dijet§ and multi-jet c3 cross-sections. 

Inclusive one-jet and dijet cross-sections have also been measured in 77 scattering at 
an e + e~ centre-of-mass energy of y/s = 58 GeV at TRISTAN c3c3 and at y/s = 130 — 



172 GeV by OPAL liSEi The E^ distribution for dijet events in the range |?? jet | < 2 
measured by OPAlEI at yfs ce = 161 — 172 GeV is shown in Fig. 12a. The measurements 
are compared to a NLO calculation of the inclusive dijet cross-sectioned which uses the NLO 
GRV parametrisation for the photon 0. The direct, single- and double-resolved parts and 
their sum are shown separately. The data points are in good agreement with the calculation 
except in the first bin where theoretical and experimental uncertainties are large. 

To study the sensitivity to the choice of parametrisation for the parton distributions 
of the photon, OPAL has also measured the inclusive dijet cross-section as a function of 
\rp ct \ for events with a large double-resolved contribution obtained by requiring < 0.8 
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Figure 12: a) The inclusive e + e~ dijet cross-section as a lunction of for jets with l^ '! < 2 using a 
cone size R = 1. b) The inclusive dijet cross-section as a function of \rf e | for jets in mainly double-resolved 
event with E 1 ^ > 3 GeV using a cone size R = 1. 



(Fig. 12 b). The variables ar 1 for the two incoming photons are defined in the same way 



7 

as the 7p variable x 7 . Ideally, for direct events without remnant jets xt, = 1 and x~ = 1, 
whereas for double-resolved events both values x+ and x~ are expected to be much smaller 
than 1. 

The inclusive dijet cross-section predicted by the two LO QCD models PYTHIAi and 
PHOJETe3 differ significantly even if the same photon structure function (here LO GRV) is 
used. This model dependence reduces the sensitivity to the parametrisation of the photon 
structure function. Different parametrisations were used as input to the PYTHIA simula- 
tion. The LO GRV0 and SaS-lD parametrisations describe the data equally well, but 
LAC1 e3 which increases the cross-section for gluon-initiated processes overestimates the 
inclusive dijet cross-section significantly. As in the case of 7p scattering a correct treatment 
of multiple interactions is important. The PYTHIA cross-sections with and without MI 
using LAC1 differ by more than a factor of two. 



3.3 Effective parton densities 

Following a procedure developed by Combridge and MaxwellS, HI has measured the parton 
distributions of the photon, /j/ 7 , using dijet events from 7p interactions c3. This method 
is based on LO matrix elements. It exploits the fact that in the studied range 0.2 < 
x 1 < 0.7 the dijet cross-section is dominated by the parton scattering processes qq' — >qq', 
qg— >qg and gg— >gg. The shapes of the angular distributions of the squared matrix elements 
\Mij(cos(9*)\ 2 is similar for all these sub-processes and the rates only differ by the ratio 
9/4 of the colour factors. The matrix element can therefore be approximated by a "Single 
Effective Subprocess" matrix element, Mses> an d the parton densities are combined into 
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Figure 13: The LO effective parton distribution a 1 x J f J (x J ,p^) 
as a function of the squared parton transverse momenta, p\. 



c 7 (/ 9 / 7 (x 7 , + (9/ 4)/ g/7 (x 7 , Pt)) 



effective parton density functions: 

/ 7 (z 7 ,Pt) = (/^(^iPt) + /g/ 7 (x 7 ,PT)) + ^fg/fi^iPr) 
n f 

fp(x 7 ,PT) = {fq/p( x V&) + /?/p( x 7>Pt)) + jfg/p(ph^) 

where the sum runs over the number of quark flavours, nf, and is the transverse mo- 
mentum of the final state parton. The dijet cross-section is then replaced by 

- - OC Y ^^.PTl^fe^T) |M SE s(cOSr)| 2 . 

dxjdxpd cos 0* ^-rr* x 7 x p 

The effective parton distribution / 7 is extracted from the data by measuring the double- 
differential dijet cross-section d 2 a/dx 7 dlogpx- 

The HI measurements of the effective parton distributions in two different x^ ranges 
(0.2 < Xj < 0.4 and 0.4 < x 7 < 0.7) are shown Fig. as a function of p\. The \ogp\ 
dependence of / 7 can now be compared to the QCD evolution using the inhomogeneous 
DGLAP equations. This is done using the LO GRV parametrisationst3 for the pion and the 
photon. Assuming that the parton distribution of the pion and the p are similar, the purely 
hadronic (VDM) part of the parton distribution functions is estimated by scaling down the 
pion distribution with a VDM factor related to the photon-p conversion probability. The 
VDM picture fails to describe the data. 

The LO GRV photon parametrisation is constructed to be purely hadronic at very 
low scales. Due to the pointlike 7 — > qq term in the DGLAP evolution used to evolve the 
photon distributions to large scales, one obtains the logarithmic rise observed for the photon 
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structure function K7(x,Q 2 ) in ej scattering (see Sect. and Fig. ||b). HI compares the 
LO GRV distribution of the photon for all partons and for quarks only to their data and 
both the rise and the normalisation of the data are well reproduced. At lower x 7 the 
gluon distribution contributes about half, whereas in the range 0.4 < x 7 < 0.7 the quark 
distribution dominates as expected. 



3.4 Leading order gluon densities 



The gluon density in the photon has been extracted directly from data by HI using a two- 
step procedure. First the distribution of the momentum fraction x 7 is unfolded from the 
data. In the second step the x 7 distribution of the gluon initiated processes is obtained by 
applying correction factors which are based on the process definitions used in the Monte 
Carlo. Since these process definitions are valid in LO, only LO gluon densities can be 
extracted. 

HI has originally performed this measurement using dijet eventsE3. This analysis suffers 
from two large systematic uncertainties related to the jet finding: the knowledge of the 
energy scale of the HI calorimeter and the effect of multiple interactions. In a new analysis 
of the gluon density, HI has avoided these uncertainties by using events with high transverse 
momentum hadrons instead of dijet events @. 
The price to be paid is a stronger sensitivity o ; 
to hadronisation effects. Events are required to ^ 
contain at least one charged particle with trans- g> 
verse momentum pr > 2.6 GeV/c and pseu- 
dorapidity \n\ < 1. The variable x™ c is re- 
constructed by summing over all hadrons with 
p T > 2.0 GeV/c and \rj\ < 1: 



^7 n 



with being the photon energy reconstructed 
from the low Q 2 tagged electron (Q 2 < 
10~ 2 GeV 2 ). 

After unfolding the LO gluon distribution 
x 7 g(x 7 )/a is shown in Fig. 14. Both analyses 
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steeply as predicted 
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Figure 14: LO gluon distribution x 1 g(x~ l ) / a of 
the photon. 



rises at low x 7 , but not as 



by the LAC1 distribution bj. ine average 
verse momentum of the final state partons. 
38 (GeV/c) 2 for the hadron analysis. 



is given by the average squared trans- 
It is 75 (GeV/c) 2 for the dijet analysis and 



3. 5 Virtual photon structure 

The measurement of the photon structure function of virtual photons in e + e~ collisions 
requires the detection of both scattered electrons. Only PLUTO has published such a 
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measurement for a virtuality Q 2 of the probing photon of 5 GeV 2 and a virtuality P 2 of the 
target photon of 0.35 GeV 2 . More double-tag measurements are to be expected soon from 
the LEP2 data. 

HI has studied the structure of the virtual photon by measuring the Q 2 dependence of 
jet production E>3. This introduces a new scale in addition to the transverse energy of the 
jet which is the relevant scale for real photoproduction (Q 2 = 0). In the kinematic range 
Q 2 » (E?p) 2 the "classical" picture of deep-inelastic scattering (DIS) can be applied where 
the pointlike virtual photon probes the parton content of the proton (Fig. p~5|a) . The jet 
transverse energy Ej, is here calculated in the 7*p centre-of-mass system. However, a small 
fraction of the events contains jets with (E^) 2 > Q 2 . In this case the process can be viewed 
as probing the structure of the virtual photon using the partonic structure of the proton. 



photon remnant 



(a)' 



jets 



proton remnant 




2 *2 
Q »E t 

proton is probed 

(deep-inelastic) 

"direct" 




H1 Preliminary 



photon is probed 
"resolved" 



(b) 



5 < E, < 7 GeV 




Figure 15: a) Schematic drawing of the two kinematic regimes studied in this analysis, b) Inclusive 7*p jet 
cross-section as a function of Q 2 for jets with 5 < E^ < 7 GeV and — 2.5 < if < —0.5 (in the 7*p frame). 



HI has measured the inclusive jet cross-section <t 7 * p for the process 7*p — > jet + X as a 
function of Q 2 for different Ej, ranges. In Fig. p)b the cross-section cr 7 * p (Q 2 ) is shown for 
jets with 5 < Ej, < 7 GeV in comparison to various MC models. LEPTO just simulates 
deep-inelastic scattering, i.e. it contains no resolved photon processes. RAPGAP contains 
both deep-inelastic scattering and resolved processes. Within RAPGAP the SaS-2D parton 
distributions for the photon is used which include a model for the Q 2 suppression for both 
the non-perturbative VDM and the perturbative anomalous part of the virtual photon. 

The DIS picture completely fails to describe the data for scales Q 2 « (E^,) 2 , but it 
approaches the data for Q 2 (E£f « 36 GeV 2 , The RAPGAP model with SaS-2D parton 
distributions is in good agreement with the data over the whole Q 2 range. This observation 
is complemented by a measurement of the fraction of the photon's energy which is assigned 
to the photon remnant. This fraction decreases with increasing Q 2 as expected if the 
resolved photon component is suppressed with increasing Q 2 . 
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4 Prompt photon production 

Prompt photon production in 7p interactions at HERA is dominated by the LO direct 
Compton process 7q — > 7q and the resolved processes qg — ► q7 and qq — > gy. "Prompt" 
means that these photons are not produced in the fragmentation process or by particle 
decays. The resolved process could be used to constrain the quark content of the photon at 
medium x 7 , complementing the -FrJ measure- 
ments from e7 scattering. The advantage of 
having a clean final state without hadronisation 
uncertainties is however largely compensated by 
the small cross-section. 

ZEUSE3 has measured the production cross- 
section for isolated photons with a transverse 
energy of 5 < E^ < 10 GeV and in the pseu- 
dorapidity range —0.7 < rf 1 < 0.8, in associa- 
tion with a jet of Eji > 5 GeV in the range 
— 1.5 < r/ jet < 1.8. Within a total error of more 
than 25 % the cross-section was found to be in 
good agreement with NLO calculations by Gor- 
donE3 using the GS and NLO GRV parton den- 
sities for the photon. ZEUS has measured x 7 by 
summing over the energies E and longitudinal momentum components p z of the photon and 
the calorimeter cells which are part of the jet 

meas _ S 7 Jet(-^' ~ Pz) 

7 2y JB E e ■ 

The variable ?/jb = J2h(E — Pz)/2E e is calculated using the calorimeter cells associated to 
the final-state hadrons. The x™ eas distribution of the ZEUS prompt photon signal is shown 
in Fig. [if]. A clear peak at x™ cas > 0.8 is observed. Using the LO Monte Carlo generator 
PYTHIA about 75 % of the events in this region can be attributed to direct Compton 
processes with x 7 = 1 with the remaining 12 % due to resolved events and 13 % due to 
radiative events where one of the outgoing quarks in a dijet event radiates a photon. More 
data are needed to obtain a more quantitative constraint for the quark distribution in the 
photon. 

5 D*^~ production in 77 and 7p interactions 

In a similar manner to jet production, open charm production in 77 and 7p collisions can 
also be used to constrain the parton content of the photon. The charm .production cross- 
sections have been calculated in NLO for 77 @'H ancWp interactions EjHH. The NLO 
calculations are either done in the so-called massive E3£3 or in the so-called massless P 6 lr 8 lr 9 I 
scheme. 

In the massive scheme the mass m c of the charm quark sets the scale for the pertur- 
bative QCD calculation. The cross-section is factorized into the matrix elements for the 
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Figure 16: x™ cas distribution of prompt photons 
measured by ZEUS. 



FREIBURG-EHEP-97-19 



18 




Figure 17: a) Cross-section for the process e + e — > e + e ccX as a function of the electron beam energy; b) 
cross-section for the process 7p — > ccX as a function of the 7p centre-of-mass energy W lv ). 



production of heavy quarks and the parton densities for light quarks (q) and gluons. This 
'massive' approach is expected to be valid if the transverse momenta pt of the charm quarks 
are of the same order, p^ ~ m c . At LEP energies only the direct process 77 — > cc and the 
single-resolved process gq^ cc are important. The relevant processes at HERA are gq^ cc 
and gg^ cc. The number of gluon initiated processes depends very much on the parametri- 
sations of the parton densities used. Again, only the sum of the event classes called resolved 
and direct is well defined in NLO. In the 'massless' scheme 
charm is considered as one of the active flavours in the 
parton distributions like u,d,s. This scheme is expected to 
be valid for pr >> Tn c . 

The cleanest method to tag open charm is the re- 



> 



construction of D* 



D 7r + decays. Due to the small 



branching ratios of the D° into charged pions and kaons, 
this method is statistics limited. ALEPH has measured 
the charm cross-section o"(e + e _ — * e + e~ccX) using_33 ± 8 
D*^ mesons reconstructed in their LEP1 data with 
Prj?* > 2 GeV. L3 has measured the charm cross-section 
in 77 interactions at LEP1 and LEP2 by tagging muons 
from semi-leptonic charm decays in the momentum range 
2 < Pu< 7 GeV/c at LEP1 and 2 < Pfl < 15 GeV/c at 
LEP2 El. The efficiency to tag muons is less than 10 
leading to large systematic and statistical uncertainties. 

The cross-section for the process e + e~ — > e + e~cc as a 
function of the beam energy is shown in Fig. 17 a. The ex- 
perimental results for various charm tagging methods used 
by pre-LEP experiments have been extrapolated to obtain 
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a total charm cross-section i. The upper band shows the fnll NLO charm cross-section 
calculated in the massive scheme by Drees et alH and the lower band the contribution from 
the Born term direct process (Quark Parton Model, QPM). The upper edge of the band 
is obtained by setting m c = 1.3 GeV with a scale \i = m c and the lower edge by setting 
m c = 1.7 GeV with fi = \^2m c . The data points obtained from the TRISTAN and JADE 
measurements cluster around the higher edge of the the massive NLO calculation which 
uses the GRV parametrisation. 

HlfH and ZEUS0H have measured D* ± production cross-sections in 7p interactions. 
Both experiments have derived a cross-section for the process 7p — > ccX as a function of 
the 7p centre-of-mass energy W 7P BEl In Fig. 0b this cross-section is compared to the 
results of lower energy experiments and to a massive NLO calculation using the NLO 
GRV parametrisation. Within the large errors the massive NLO calculations using the 
GRV parton distributions are in good agreement with the total charm production cross 
sections measured at LEP and HERA. The cross-section depends very much on the gluon 
distribution in the photon, e.g. the LAC1 parametrisation gives a much larger cross-section. 

The extrapolation of the D*^ production cross-section to a total cc cross-section has 
large theoretical and experimental uncertainties. These are avoided if differential distribu- 
tions are measured. ZEUsEEl has measured the p®* distribution (Fig. [l8| ) and compared 
it to a massive NLO calculation 153 using the GRV parametrisation for the photon and to 
a massless calculation E2. The normalisation of the massless calculation which uses similar 
parameters as in the case of the massive calculation with a charm mass m c = 1.5 GeV is in 
better agreement with the data. 



6 Total cross sections 

The total cross-sections a for hadron-hadron and 7p collisions are well described by a Regge 
parametrisation of the form a = Xs e + Ys~ v , where yfs is the centre-of-mass energy of the 
hadron-hadron or 7p interaction. The first term in the equation is due to Pomeron exchange 
and the second term is due to Reggeon exchange H. Assuming factorisation of the Pomeron 
term X, the total hadronic 77 cross-section <r 77 can be related to the pp (or pp) and 7p 
total cross-sections at high centre-of-mass energies W = -^s where the Pomeron trajectory 
should dominate: 

a 77 = (4) 
o-pp 

A slow rise of the total cross-sections with energy is predicted, corresponding to e ~ 0.08. 
This rise can also be attributed to an increasing cross-section for parton interactions leading 
to mini-jets in the final statetiZI. 

Before LEP ct 77 (IU) has been measured by PLUTO 0, TPC/27H and MD10. These 
experiments have measured at 77 centre-of-mass energies W below 10 GeV before the 
onset of the high energy rise of the total cross-section. Using LEP data taken at y^ee = 
130 - 161 GeV L30 has demonstrated that a, w(W) is consistent with the universal Regge 
behaviour of total cross-sections in the range 5 < W < 75 GeV which was also observed in 



7P scattering at HERA. The L3 measurement is shown in Fig. 19 a together with an OPAL 
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Figure 19: Total cross-section ol the process (a) 77 — > hadrons as a function of W = W 7 
the process (b) 7p — > hadrons as a function of W lp = \/s . 
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measurement @ in the range 10 < W < 110 GeV using data taken at \/s ee = 161 — 172 GeV. 
The observed energy dependence of the cross-section is similar, but the values for cr 77 
are about 20 % higher. The errors are strongly correlated between the W bins in both 
experiments. Furthermore, L3 has used the Monte Carlo generator PHOJET e3 for the 
unfolding, whereas OPAL has averaged the unfolding results of PHOJET and PYTHIA. The 
unfolded cross-section using PHOJET is about 5 % lower than the central value. In both 
experiments the cross-sections obtained using PHOJET are lower than the cross-section 
obtained with PYTHIA. The origin of the remainingdiscrepancy is not yet understood. 

Based on the Donnachie-Landshoff (DL) model tZI, the assumption of a universal high 
energy behaviour of 77, 7p and pp cross-sections is tested. The parameters X and Y are 
fitted to the total 77, 7p and pp cross-sections in order to predict <7 77 via Eq. |||. This is 
done assuming that the cross-sections can be related at \/s = v^-yp = V^pp • The process 



7P 



PP 



dependent fit values for X and Y are taken from Ref. 73 together with the values of the 
universal parameters e = 0.0790 ± 0.0011 and r\ = 0.4678 ± 0.0059. This simple ansatz gives 
a reasonable description of the total 77 cross-section o" 77 . Schuler and Sjostrand give a 
total cross-section for the sum of all possible event classes in their model of 77 scattering 
where the photon has a direct, an anomalous and a VDM component. They consider the 
spread between this prediction and the simple factorisation ansatz as conservative estimate 
of the theoretical band of uncertainty. The prediction of Engel and Ranft e3 is also plotted 
which is implemented in PHOJET. It is in good agreement with the L3 measurement and 
significantly lower than the OPAL measurement. The steeper rise predicted by Engel and 
Ranft is in agreement with both measurements. 



7 Photon fragmentation function 

Closely related to the parton distributions of the photon are the fragmentation functions 
D^ s (z, M 2 ) of quarks and gluons into photons with z being the fractional momentum 
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carried by the photon. The photon fragmentation functions (z, M 2 ) are measured in 
e + e~ annihilation at LEP. The time-like scale M 2 is therefore given by the e + e~ centre-of- 
mass energy. 

As in the case of the photon structure function, the photon fragmentation function 
is fully calculable in perturbative QCD for asymptotically large M 2 due to the pointlike 
coupling of the photon to qq pairs. At experimentally accessible values of M 2 the non- 
perturbative contributions are still large and have to be taken into account. Bourhis et al.El 
have calculated new fragmentation functions with a full treatment of the Beyond Leading 
Logarithm (BLL) corrections to the perturbative part and with a VDM input for the non- 
perturbative part. 

A first measurement of the quark-to-photon 
fragmentation function has been published by 
ALEPH using isolated photons as well as photons 
which were reconstructed inside hadronic jets 0. It 
was pointed out by the authors of Ref. 75 that these 
measurements cannot be compared directly to the 
fully inclusive calculations of the photon fragmen- 
tation functions due to phase space restrictions im- 
posed by the jet algorithm used. OPALO has there- 
fore measured the fully inclusive energy spectrum 
for photons with energy E~ > 10 GeV in hadronic 
Z° decays. The huge background from photons due 
to hadron decays (e.g. 7r° — > 77) is subtracted us- 
ing a fitting method. The data shown in Fig. |2^ are 
in agreement with the models by Bourhis et al El, 
Duke and OwensEl and by Gliick et al.El. The dif- 
ferences between the GRV LO calculation and the 
higher order (HO) calculation with and without the 
non-perturbative corrections (HOpl) is much smaller than the experimental errors. 

8 Conclusions 

Interactions of photons via quantum fluctuations can be described using a structure func- 
tions formalism. The hadronic structure function F^x^Q 2 ) of the photon is measured in 
e7 scattering at LEP in the range x > 10 -3 and 1 < Q 2 < 10 3 GeV 2 . The logarithmic rise 
of F2 with Q 2 for medium x and large Q 2 is observed as predicted by perturbative QCD. 
At low x LEP will be able to study the region where the onset of the rise of F2 is expected 
from the HERA data on the proton structure function. 

The parton content of the photon is also measured in 77 and 7p interactions at LEP 
and HERA. Jet production and high px hadron production are especially sensitive to the 
gluon distribution in the photon. 

The GRV, SaS, GS and LAC parton distributions of the photons are currently the 
most widely used parametrisations. Parametrisations with a large gluon distribution like 
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LAC1 are disfavoured by the data. The other distributions are consistent with most of the 
available measurements. 

HI has presented new results on the structure of virtual photons from the Q 2 dependence 
of jet production. More information about the interactions of photons at low and medium 
Q 2 are to be expected in the near future from the measurement of double-tagged 7*7* 
events at LEP. 

First measurements by L3 and OPAL of the energy dependence of the total 77 cross- 
section for hadron production show the rise characteristic of hadronic interactions. 
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